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ABSTRACT 

We have studied the determination of the running 6-quark mass, m^Mz), using Z° 
decays into 3 or more hadronic jets. We calculated the ratio of > 3-jet fractions 
in e + e~ — > bb vs. e + e~ — > q{qi (qi = u or d or s) events at next-to-leading order in 
perturbative QCD using six different infra-red- and collinear-safe jet-finding algorithms. 
We compared with corresponding measurements from the SLD Collaboration and found 
a significant algorithm-dependence of the fitted rrib{Mz) value. Our best estimate, 
taking correlations into account, is mj,{Mz) = 2.56 ± 0.27 (stat.) tgjg (syst.) ti.Jg 
(theor.) GeV/c 2 . 
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1. Introduction 



Three-jet events of the type e + e~ — > qqg provide an ideal laboratory for making pre- 
cise tests of Quantum Chromodynamics (QCD) [p]]. Since the initial state is free of 
strongly-interacting particles the experimental environment is intrinsically 'clean', and 
the process is more amenable to calculation using perturbation theory than, for ex- 
ample, multijet final states in hadron-hadron or lepton-hadron collisions. A number 
of perturbative QCD (pQCD) predictions for 3-jet dominated hadronic event-shape 
observables, for massless quarks, complete at next-to-leading order (NLO) are avail- 
able § 

One would expect the emission of gluon radiation in events containing massive 
quarks, e + e~ — > QQ (Q = b or c), to be modified relative to that in e + e~ — > qiqi (qi 
= u or d or s) events due to the restriction in phase space imposed by the non-zero 
quark mass. One would also expect such a modification to depend on the choice of 
the event-shape observable, and to be potentially relatively large for those observables 
in which the quark-jet mass enters kinematically into the definition. Recently NLO 
calculations of e + e~ — > QQg have been performed in which quark mass effects have 
been taken into account explicitly || [II], |TTJ . ^From these calculations one expects the 
size of the c-mass effects in ccg events at yfs = Mz to be rather small, at the level 
of 1%. However, for bbg events the relative size of the 6-mass effects on event-shape 
observables can be much larger, up to around 5%. Such a large effect needs to be 
taken into account in precise studies of bbg events where the experimental errors can 
be comparable with, or smaller than, this size. 

For example, tests of the flavour-independence of strong interactions involve mea- 
surements of the ratios r®(X) = X®/X of a 3-jet observable X in QQg versus qiqig 
events. Currently the experimental errors on r b (X) are of the order of 1-2%, and b- 
mass effects are clearly visible in the data |T3|. By contrast, the errors on r c (X) 
are much larger than 1% and any c-mass effects are not discernible. Hence, in recent 
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measurements the NLO massive calculations have been employed to correct r b (X) for 
the b- mass effects, so as to determine the ratio of strong couplings, a b s /a^ ds [13|, [y], [15| , 
and test the ansatz of flavour-independence of strong interactions. 

An alternative, and a priori equally valid, approach is to assume that strong inter- 
actions are flavour independent, and use the sensitivity of event-shape observables to 
mass effects to determine the 6-mass itself. In the theoretical prediction one has the 
freedom to choose the renormalization scheme which defines the quark mass. For ex- 
ample, one can write the NLO result in terms of either the perturbative pole mass Mj, 
or the 'running' mass m&(/i). The latter is defined by the modified minimal subtraction 
(MS) scheme ]16[ employed to renormalize the mass at a scale fi. At the Z° scale, 
Mz, the running mass is preferable because large logarithms of the form ln(M^/Mf) 
are absorbed in m^iMz)-, and the perturbative expansion is thus improved. The DEL- 
PHI Collaboration has recently studied the 3-jet-rate i?3(y c ), where R3 was determined 
using the Durham (D) jet-finding algorithm [17], and y c is the scaled-invariant-mass 



criterion which determines the jet multiplicity. From their measurement of r (R3) at 



y c = 0.02 DELPHI obtained 13 



m b (M z ) = 2.67 ± 0.25(exp.) ± 0.34(frag.) ± 0.27(theor.) GeV/c 2 . (1) 

Such 3-jet-event observables have been used for many years to determine a s (Mz) 
from inclusive-flavour e + e~ annihilation events [|18j. Though the a s value obtained by 



fitting a NLO pQCD calculation to any one measured observable can have quite a small 
experimental error, < 0.001 for some observables, there is a strong dependence of the 



fitted a s (Mz) value on the choice of observable fl9|, pPf. This spread of values leads to a 
large and dominant uncertainty on a s [Mz) determined with this technique [18] . Since 
non-perturbative effects are supposedly taken into account in these measurements, usu- 
ally by applying corrections based on well-tested hadronisation models [^TJ, a consistent 
description within the framework of pQCD is viable only if one postulates large (and 
uncalculated) next-to-next-to-leading-order (NNLO) contributions to the observables. 
Hence, in this picture, the spread in a s [Mz) values determined at NLO results from 



the omission of the uncalculated higher-order terms. Furthermore, it can be argued 
that a strong dependence of a NLO calculation on the renormalisation scale is generally 
a sign of large NNLO contributions. Such a dependence is indeed observed for most of 
the observables |T^, |2^, |22|] , and supports the previous interpretation, though there is 
little consensus on a procedure for quantifying the scale-dependence of measurements 
of a s (M z ). 

The DELPHI determination of rrib(Mz) (Eq. |IJ) is based on a ratio of 3-jet-event 
observables calculated at NLO. Given this apparently very precise result derived from 
one observable, it is interesting to consider the possible effect of NNLO contributions. 
Naively one might expect any potentially sizeable effects in the numerator and denom- 
inator largely to cancel. However, a residual uncertainty at only the 2% level on r b 
corresponds (Section 3) to a 0.5 GeV/c 2 error on m^Mz), which is comparable with 
the quoted total error on the DELPHI measurement. For the purpose of investigation 
we have studied the extraction of rrib(Mz) from the ratios r b {R^), where R3 was deter- 
mined using six different infra-red- and collinear-safe jet-finding algorithms. As in the 
case of a s (Mz) measurements using such observables, the study of an ensemble of re- 
sults from different observables, all calculated at NLO, may uncover systematic effects 
relating to the uncalculated NNLO contributions. We used the Durham and Geneva 
(G) p3| schemes, and the E, E0, P and P0 variations of the JADE algorithm [23] to 
evaluate the 6-mass- dependent NLO pQCD predictions, and compared them with the 



corresponding experimental measurements published by the SLD Collaboration [14 



In Section 2 we outline the theoretical framework and briefly describe the NLO cal- 
culations used here. In Section 3 we compare the calculations with the data and extract 
values of m&(Mz) using each jet algorithm in turn. We compare the results obtained 
using the different jet algorithms, and discuss the systematic uncertainties. In Section 
4 we combine these results to obtain our best estimate of the central rrib(Mz) value 
and error by taking correlations into account. These results supercede the preliminary 



results presented in W2 
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2. Theoretical Framework 



In this section we describe the computation of the double ratio: 

r b (y c ) = R b (y c )/Rf s (y c ), (2) 

where we define R\{y c ) to be the fraction of events classified as containing 3 or more 
jets with a particular jet-finding algorithm at a particular y c value. The event flavour 
is defined by the flavour of the primary quarks that couple to the Z°. This definition 
means that events of the type Z° — > qiqig — ^ qiqibb are classified as light-quark events 
[] R b (y c ) and R^ ds (y c ) can be written to NLO accuracy: 

RliVc) = ^A«(y c ) + (B"(y c ) + C*(y c )) + O (afj , (3) 

where the coefficients A (B) represent the LO (NLO) contribution to 3-jet production, 
and the coefficients C represent the LO contribution to 4-jet production. Thus we have 
(we suppress here the argument y c ): 



A b a s (B b + C b B uds + C uds A b 

+ 7T- 



) + O (« s 2 ) . (4) 



j^uds 27]- \ j^uds j^uds j^uds J 

Some comments are in order about the generalization of observables originally de- 
fined for massless quarks to account for quark masses. Because of the fact that, for 
massless quarks, energy and the modulus of the three-momentum can be interchanged, 
the generalization to massive quarks contains a certain degree of freedom. For massive 
quarks one should use a definition which is infrared safe and does not distinguish be- 
tween a hard parton and the hard parton substituted by two collinear ones, or a hard 
parton and soft parton. (See for example |25| where this is discussed for the example 



of the thrust.) While this is not a problem for the Durham, Geneva, and E algorithms 
1 This distinction is possible only because the contribution to i?| from the interference of the 
amplitudes for Z a — > qiqig — * qiqibb and Z° — * bbg* —> bbqiqi is negligible. In fact the contribution 
cancels exactly in i?g if one neglects the masses of the u,d,s,c quarks and the sum over the light 
quarks qi includes only full isospin doublets. 
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one should be careful when using algorithms which use a recombination procedure that 
keeps the recombined four-momentum massless, in particular the EO, P, and PO algo- 
rithms. Replacing a massive quark by a hard parton and a soft one, or two collinear 
partons, would yield a massless four-momentum vector after the recombination, in- 
stead of getting the original four-momentum back. These algorithms hence 'feel' soft 
and collinear partons if they are naively applied for massive quarks. (Fortunately this 
problem does not occur if further splitting is considered.) 

One possibility would be to modify these algorithms in some fashion so as to keep the 
hadronisation corrections small; this was in fact the original motivation for developing 
the Durham and Geneva algorithms. Here we preserve the definition of the algorithms 
so as to be able to compare with the published results of the SLD analysis. The SLD 
measurements are corrected to the parton-level, where the naive jet definition is applied 
even if one of the four-momenta is massive. We are hence able to extract consistently 
the b mass from these measurements obtained with the E0,P, and PO jet algorithms by 
treating the b quark as massive in our theoretical prediction. 

A. 3-jet Contributions 

For 3-jet production the LO massless (A uds ) and massive (A b ) |26| as well as the NLO 
massless (B uds ) || coefficients are well known. In order to calculate the massive NLO 
coefficients B b we need the matrix elements for 

e+e-^(j\Z*)^bbg (5) 

to order a 2 s , as well as the matrix elements for the parton processes 

e + e" -> (7*, Z*) -> bbgg, bbq { q h bbbb. (6) 

In the calculation of the virtual corrections to Eq. |5| both ultra-violet (UV) and infra- 
red (IR) singularities are encountered. The UV singularities are removed by the usual 
renormalization procedure of the mass parameter and the QCD coupling a s . 
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The IR singularities are cancelled by the real contributions from the processes listed 
in Eq. @. It is worthwhile adding some remarks about this cancellation. Today it is 
more or less standard to regulate the IR divergences in the framework of dimensional 
regularization. To cancel the divergences in the virtual corrections one must then 
integrate the real contributions over some regions of phase-space in d dimensions. More 
precisely one must integrate over the regions where a gluon is soft or two massless 
partons are collinear. In general this would be a formidable task. Therefore several 
techniques have been developed in the past (see for example ||, |^, ^ [27]]) to simplify 
this problem using the general factorization properties of QCD amplitudes. In the 
calculation reported in |J on which the current paper is based the so called phase- 
space-slicing method |3| was used. The same is true for the results presented in |TI| on 
which the DELPHI analysis |13| is based. In the calculation given in Jll| an alternative, 
the so-called subtraction method, was used. 

In the simplest version of the phase-space slicing method one separates the 'soft' 
and 'collinear' regions (often called 'unresolved regions') from the rest of the phase- 
space ('resolved regions') by demanding a minimal invariant mass-squared s m i n for all 
pairs of partons. In the soft and collinear regions the squared matrix elements can be 
approximated by the use of the soft and collinear factorization which is valid in the 
appropriate limits. After this simplification the relevant part of the squared matrix 
elements can be integrated analytically in d dimensions. The phase-space integration 
over the resolved regions can be done numerically in four dimensions. In the case of 
massive quarks the phase-space slicing method must be modified, although the basic 
features are the same. In particular, the 'slicing' between the soft/collinear regions 
and the regions where all partons are hard can still be parametrized in terms of one 
variable s mm . 

The approximation used in the unresolved region is only valid for small values of 
Smin- On the other hand for small s m ; n large cancellations between the numerically 
integrated and the analytically integrated parts will arise leading to possible errors in 
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the sum of the two. Note that the artificial cut parametrized by s min is not related 
to any physical cut. Thus the theoretical prediction must be independent of s min . In 
practice, for the NLO coefficient B b this will be true up to corrections of the order of 
s min /s. With the value s min = 0.5 GeV 2 , which we have used in our calculation, the 
systematic error in B b due to the phase-space slicing method is negligible compared 
with the numerical error due to the numerical integration, which is itself negligibly 
small. 

Note that although free of collinear singularities in the case of massive quarks, one 
must also include the contribution of the reaction e + e~ — > (7*, Z*) — > bbbb to the 3-jet 
rate. The collinear singularity appearing for massless quarks shows up as \n{m b ) for 
massive quarks. Only if this contribution is included is the 3-jet contribution free of 
collinear singularities in the limit of vanishing 6-quark mass. The logarithm cancels 
against the logarithm coming from the virtual corrections to the gluon propagator due 
to a 6-quark loop. Needless to say, this reaction contributes also to the 4-jet rate. 

Note also that in the reaction e + e~ — > (7*, Z*) — > bbgg collinear singularities (due 
to collinear configurations of a (anti-)quark gluon pair) that appear for massless quarks 
give rise only to logarithms of the mass. This explains why it is difficult to calculate 
i?3 numerically for small values of m b without further analytical work. In the virtual 
correction these logarithms are explicit. In the real contribution they must be obtained 
numerically from the phase space integration, which becomes more difficult for smaller 
masses due to numerical instabilities. 

For technical reasons it is easier to perform the calculation of r b first in the pole 
mass scheme, and switch to the running mass afterwards. The relation between the 
pole mass M b and the MS mass m b (fi) reads 



M b = m b (fx) + Sm b (ii), 



(7) 



where, to order a, 
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This implies the following relation between r in both mass renormalization schemes 28 



The mass dependence of A b can be written as A b (m b ) = A uds + 5A b (m b )m b 1 / s, which 
defines the function 5A b (rrib). For m 2 . <C s, 5A b depends only weakly on m b . Ignoring 
this residual mass dependence we have 

1 (L4 6 (m m (//)) 1 

-5m b {n)- 



A uas--°^> dmb{fl) - A u ds ^ [y/mlifi) + 2m b (ri5m b (v)) - A b {m b ^))\ , (10) 

and we use the r.h.s. of ([II]) to convert the results for r b from the pole mass to the MS 
mass renormalization scheme. This excellent approximation avoids the calculation of 
the derivative dA b /dm b for each algorithm. 



B. 4-jet Contributions 

Both the massless [[J and massive |29j LO 4-jet fraction contributions (C) are well 
known. Recently, the massless 4-jet fraction has been computed to NLO f30| . These 
corrections, which are of order and therefore not included in our prediction for r b , 
can change, depending on the jet algorithm, the values of the massless C coefficients by 
up to 100%. Note that part of these large NNLO corrections to r b will cancel between 
the massless and (yet unknown) massive O (af) C coefficients entering (|]). 



C. Calculation of r h 

The calculation of r b was performed, for each of the six jet algorithms at the optimal y c 
values (discussed below), for m b (Mz) values in the range 2.0 < m b (Mz) < 5.0 GeV/c 2 . 
These predictions are shown as points in Fig. |l|. For the P and P0 algorithms we also 
calculated a point at m b (Mz) = 1 GeV/c 2 so as to constrain better the extrapolation 
r b — y 1 as m b (M z ) — > 0. The renormalization scale /i was set to ^fs and we used 
a s (Mz) = 0.118. The dependence of r b on the renormalization scale fi is trivial if r b 
is expressed in terms of the pole mass M b ; it enters only through the running of a s . 

9 



An additional fi dependence is introduced if one switches to the running mass rrib(fi) 
by using Eq. (H). The theoretical uncertainty on m^Mz) due to the choice of the 
renormalization scale will be discussed in section 3. 
The function 



„. , m 2 „m 2 . ( m 2 \ m 4 

f(m) = 1 + + In ( — ] +7— 



, , (11) 

S S \ S J 

where a, f3 and 7 are free parameters, was fitted to these points. The ansatz (|TT|) can be 
jusified as follows: 1) As m -> 0, the massive fraction R\ approaches the massless oneQ 
R^ ds ; 2) since m 2 <C s, it is a very good approximation to keep only the leading terms 
in m 2 j s. The fitted parameter values are listed in Table |1J the functions are shown in 
Fig. [I] and provide a good description of the mass dependence of the calculations. 



Algorithm 


Vc 


a 


P 


7 


E 


0.040 


207.6 


16.10 


-13029.9 


E0 


0.020 


42.2 


-3.58 


-3881.3 


P 


0.020 


194.6 


26.17 


-3855.9 


P0 


0.015 


206.0 


26.73 


-1317.5 


D 


0.010 


79.3 


17.16 


-4610.8 


G 


0.080 


-89.6 


-11.04 


3229.9 



Table 1: Fitted parameters of the function /(to) = 1 + am 2 / s + (3{m 2 / s) ln(m 2 /s) + 
7m 4 / s 2 for each jet algorithm. 



It can be seen that the m6(Mz)-dependence varies according to the jet algorithm. 
For m b (Mz) > 2.0 GeV/c 2 , r b > 1 and the slope is positive for the E, E0, P and P0 
cases, whereas r b < 1 and the slope is negative for the D and G cases. This can be 
understood qualitatively in terms of two competing physical origins. First, the non- 
zero 6-mass tends to cause a phase-space suppression of gluon emission relative to the 
2 This is true up to differences induced by triangle diagrams pi) , which lead to deviations of the 
r b {m = 0) from 1 of less than 0.1%. 
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massless quark case, implying r b < 1. Second, for a given kinematic configuration, the 
large 6-mass tends to enhance the invariant mass of a local quark-gluon pair relative 
to the massless quark case. Since the JADE family of jet algorithms is based on a 
clustering metric that is closely related to invariant mass, for fixed y c the two partons 
are more likely to be resolved as separate jets when the quark is massive, implying 
r b > 1. By contrast, the clustering metric used in the Durham and Geneva algorithms 
is less sensitive to this kinematic effect, the phase-space suppression dominates, and 
r b < 1. For increasing values of y c one expects both effects to diminish in importance 
and r b — > 1. For the D algorithm this has been observed in the DELPHI study fl"3| . 



3. Extraction of the 6-Quark Mass 



We used measurements of r published |L4| by the SLD Collaboration. These mea- 



surements are based on a sample of 150,000 hadronic Z° decays recorded between 
1993 and 1995, for which the original 120-million pixel CCD vertex detector was used 
for event flavour separation. SLD has subsequently recorded a further 400,000 Z° 
decays with a new 30 7- million-pixel vertex detector, and it would be straightforward 
to repeat the present analysis when the new data are made available. Though not as 
statistically powerful as the DELPHI result for the D jet algorithm, the SLD published 
data include results for the six different jet algorithms D, G, E, E0, P and P0, and are 
hence suitable for this study of possible observable-dependent systematic effects. 

Full details of the experimental procedure are given in [|14]]. Briefly, e + e~ — > hadrons 
events were selected, and a flavour-tagging algorithm was applied to select samples of 
events of primary b, c, and uds quark flavour. The algorithm was based on the mass 
and momentum of secondary decay vertices reconstructed using the vertex detector. 
Light-quark (uds) events rarely contain reconstructed secondary decay vertices, and 
these typically result from strange particle decays and are of low mass. Conversely, 
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bb events typically contain high-mass vertices from 5-hadron decays. The purity of 
the fe-tagged (tids-tagged) event sample was 90% (91%) respectively. 

Each jet-finding algorithm was applied in turn to the uds- and 6-tagged samples and, 
for each algorithm, the ratios (Eq. were formed. The ratio is an attractive quantity 
as many of the experimental and theoretical systematic uncertainties effectively cancel. 
Each ratio was then corrected for the effects of detector acceptance and resolution, the 
bias of the flavour tag to select preferentially 2-jet rather than 3-jet events, the flavour 
compositions, and hadronisation effects. For each algorithm an 'optimal' y c value was 
selected so as to minimise the combined statistical and experimental systematic error. 

The measured r b values and the associated errors are listed in Table ^| |14j]. The 
central values and statistical errors are also shown in Fig. [IJ. The set of r b values 
is not consistent with unity, which indicates that the 6-mass effects are significant. 
Furthermore, a systematic algorithmic dependence is apparent, with r b > 1 for the 
JADE family of algorithms and r b < 1 for the D and G algorithms, in agreement with 
the expectations discussed in Section 2. 



Algorithm 


Ve 




stat. 


exp. syst. 


had. 


E 
EO 
P 
PO 
D 
G 


0.040 
0.020 
0.020 
0.015 
0.010 
0.080 


1.050 
1.054 
1.048 
1.055 
0.964 
0.995 


0.026 
0.019 
0.019 
0.017 
0.023 
0.032 


+0.038 
-0.042 

+0.030 
-0.037 

+0.027 
-0.037 

+0.028 
-0.035 

+0.038 
-0.041 

+0.035 
-0.036 


+0.011 
-0.046 

+0.007 
-0.045 

+0.002 
-0.026 

+0.007 
-0.037 

+0.001 
-0.006 

+0.020 
-0.008 



Table 2: SLD measured values and errors of r b . 



For each jet algorithm, by comparing the theoretical curve in Fig. p] with the SLD 
data, one can read off the preferred rrib(Mz) value. The central values are listed in 
Table |3|. In each case upper and lower statistical errors were evaluated from the crossing 
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points of the error band with the theoretical prediction, except in the case of the G 
algorithm, for which the upper statistical bound is consistent with nib = 0; in this case 
an error equal to the central value was assigned. Each experimental systematic error 



on r b 



14 1 was similarly transformed into a systematic error on mt,(Mz) and the sum 
in quadrature is listed in Table [3[ Hadronisation uncertainties jy] were evaluated in a 
similar fashion and are listed in Table |3|. 

Additional theoretical uncertainties were investigated by varying the value of a s (Mz) 
within the range 0.115 <u s (Mz) < 0.121. The corresponding changes in m^iMz) were 
at the level of ±(10 — 20) MeV/c 2 . The renormalisation scale was also varied within the 
range Mz/2 < /i < 2Mz- The corresponding changes in vn^iMz) were at the level of 
±200MeV/c 2 or less. For each algorithm these uncertainties were added in quadrature 
to define a theoretical uncertainty, which is listed in Table 

The six measured 6-quark masses range from 2.3 to 4.1 GeV/c 2 , with an r.m.s. 
deviation of 0.7 GeV/c 2 ; this scatter is larger than one might expect from these data 
given the strong correlations between measurements using different jet algorithms, 
suggesting some additional source of uncertainty. In order to quantify this issue we 
evaluated the statistical correlations among the r b values determined using different jet 
algorithms. We repeated the analysis on subsets of both the data and the simulated 
data and calculated the correlation coefficients empirically. The data and simulation 
gave consistent results, and the average correlation coefficients are listed in Table 
Each has a statistical uncertainty of ±0.03. The four JADE-like algorithms show strong 
correlations with each other, in the range 0.65-0.84, as might be expected. Correlations 
between other pairs of algorithms are weaker, in the range 0.41-0.65. 
We evaluated 

X 2 = E«(r? - f, i (m b (M z )))(V-%(r b 1 - fj(m b (Mz))), (12) 

where r\ (/,) are the measured (calculated) double ratios, i,j = E, E0, P, P0, D, G, 
the error matrix is defined by V^- = CijOiO^ Cij is the correlation coefficient given in 
Table |], and <Ji is the quadratic sum of the data and Monte Carlo statistical errors 
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Algorithm 


rrib(Mz) stat. exp. syst. had. theor. 


E 


2.271 +0.488 +0.734 +0.217 +0.194 
-0.629 -0.952 -1.483 -0.189 


EO 


2.642 +0.493 +0.789 +0.187 +0.213 
-0.562 -1.082 -1.637 -0.226 


P 


4.067 +0.423 +0.616 +0.048 +0.047 
-0.504 -0.981 -0.723 -0.021 


PO 


3.728 +0.307 +0.515 +0.132 +0.056 
-0.361 -0.738 -0.911 -0.043 


D 


2.509 +1.028 +1.879 +0.287 +0.170 
-1.255 -2.001 -0.049 -0.195 


G 


2.415 +2.075 +2.761 +0.691 +0.195 
-2.415 -2.415 -2.415 -0.078 



Table 3: Values of the running 6-quark mass extracted from the SLD measurement of 
r b for each of the six jet algorithms. 



Algorithm 


E 


E0 


P 


P0 


D 


G 


E 


1.00 


0.70 


0.67 


0.65 


0.61 


0.49 


E0 




1.00 


0.84 


0.82 


0.61 


0.49 


P 






1.00 


0.71 


0.65 


0.56 


P0 








1.00 


0.52 


0.41 


D 










1.00 


0.64 


G 












1.00 



Table 4: Statistical correlation coefficients between r b measurements for each pair of 
jet-finding algorithms. 
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on r\. For values of m^iMz) around 2.9 GeV/c 2 , which is the average of the results 
shown in Table H we obtained y 2 ~ 38/6, which indicates an inconsistency among 
the results from the different algorithms. We minimised x 2 with respect to variation 
of m b (M z ) and obtained x 2 = 26/5 for m b (M z ) = 0.9 ± 0.7 (stat.) GeV/c 2 , which 
is still unacceptably high. The best-fit x 2 value is insensitive to variations of the 
within their uncertainties, and to (simultaneous) systematic shifts of the measured 
within the experimental systematic errors and hadronisation and theoretical uncertain- 
ties. The experimental systematic errors, which are dominated by uncertainties in the 
flavour composition of the samples, and the hadronisation uncertainties were assumed 
to be 100% correlated among all algorithms and were omitted from the x 2 calculation; 
theoretical uncertainties were also omitted. 

We repeated this minimisation procedure and omitted in turn the measurement 
based on each of the six algorithms. In no case did we obtain a x 2 value better 
than 12, which corresponds to a confidence level of 1.7%. We then omitted pairs of 
measurements in turn. Fits with x 2 < 9, i.e. >6% confidence level, were obtained only 
for two of these 15 cases in which both the E and E0, or both the P and P0, algorithms 
were omitted; the corresponding mass values were 2.5 and 3.5 GeV/c 2 , respectively. 
To the extent that the hadronisation and theoretical uncertainties have been properly 
estimated, we do not have a prion justification for omitting any particular algorithm(s). 
We do note, however, that algorithms in the JADE family have a significantly worse 
soft gluon behaviour than the D and G algorithms |2lJ . The former algorithms tend to 
combine soft gluons to form an 'artificial' jet at values of y c that are not small, which 
may cause large higher-order perturbative corrections even for moderate values of y c . 

The x 2 value is, however, quite sensitive to small changes in the measured r\ or 
predicted As an exercise, we have postulated an additional uncertainty of size e 
which is uncorrelated between different jet algorithms. Under the hypothesis m& = 2.94 
GeV/c 2 , which is the average of the values listed in Table |3], for e = 0.015 the x 2 value 
is 10.5; for e = 0.02 the x 2 value is 7.1, which is acceptable. For e > 0.02 a fit 
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for mb(Mz) yields a stable value of roughly 2.6 GeV/c 2 , indicating that a consistent 
rrib(Mz) value can be obtained provided that there exist additional uncertainties, un- 
corrected between jet algorithms, at the level of 2% on r b . A 2% error on the predicted 
r b corresponds to an error of ~0.5 GeV/c 2 on the extracted value of m^Mz) from a 
given algorithm, and would roughly account for the 0.7 GeV/c 2 r.m.s. deviation among 
the values in Table [| 

We suspect that the most likely source of the inconsistency among results for the 
different jet algorithms is the missing higher-order perturbative contributions to r b . As 
we have shown, these would have to be only at the level of 2% in order to resolve the 
inconsistency. Possible NNLO contributions to r b of this small magnitude are not a 
priori unexpected; 5-10% level NNLO contributions are implied by the scatter among 



a s (M z ) values determined using these and closely-related event-shape observables fll8 
which form the numerator and denominator of r b . 



4. Summary and Conclusions 

We have studied the determination of the running 6-quark mass by comparing NLO 
perturbative QCD calculations of the > 3-jet ratio r b = R\/K^ ds with data from the 
SLD Collaboration. We used six different infra-red- and collinear-safe jet-finding al- 
gorithms in order to study systematic effects. We find algorithm-dependent values of 
mb(Mz) in the range 2.3 < rribiMz) < 4.1 GeV/c 2 . The value determined using the 
Durham algorithm is consistent with that reported in [[0| . 



We quantified the statistical, experimental systematic, hadronisation and additional 
theoretical uncertainties, and attempted to obtain a best-fit m^iMz) value by minimis- 
ing x 2 , taking statistical correlations between the results from the different algorithms 
into account. We could not obtain an acceptable best-fit x 2 value unless, in two cases, 
we omitted two algorithms from consideration. In the absence of an a priori reason 
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to do this we retained all six algorithms in order to investigate possible additional 
systematic effects. We were able to obtain an acceptable value of x 2 > an d a stable 
value rrib(Mz) — 2.6 GeV/c 2 , provided that we postulated (an) additional source(s) of 
uncertainty of relative size > 2% on r b , which is uncorrelated between algorithms. We 
are unable to account for the origin of such an uncertainty, but speculate that it may 
be due to uncalculated higher-order pQCD contributions. 

We now discuss the assignment of a single value of m\,{Mz). Taking an unweighted 
average of the m b values in Table g, yields m b (M z ) = 2.94±g;|^ (stat.) t\H (syst.) t?;^ 
(had.) 1qj2 (theor.) ±0.69 (r.m.s.) GeV/c 2 , where we include the r.m.s. deviation as 
an additional error. Though well defined, this procedure does not make full use of the 
information contained in the six measurements. The \ 2 minimisation procedure does 
take into account the full statistical covariance matrix, as well as correlations in the 
systematic error and hadronisation uncertainties. Since the resulting y 2 is acceptable, 
and the fitted rrib{Mz) value is stable for any additional uncorrelated uncertainties of 
size e > 0.02, we choose e = 0.02, and obtain \ 2 — 6.5/5 with 

m b (M z ) = 2.56 ± 0.27(stat.)^;|(syst.)^g(theor.) GeV/c 2 . (13) 

We consider that this represents our best estimate of the running 6-quark mass using 
the SLD data. 

The statistical error on m^iMz) is substantially reduced by the correlations among 
the six individual r b results. The experimental systematic error is also reduced by 
the fact that a given shift of r b causes the m^iMz) values for the E, E0, P and P0 
algorithms to shift in opposite directions to those for the D and G algorithms. The 
theoretical uncertainty comprises the sum in quadrature of the hadronisation uncer- 
tainty (ii^ GeV/c 2 ), and the propagation of the uncorrelated error of ±0.02 on each 
r b (±0.42 GeV/c 2 ). Variations of a s and /i contribute an uncertainty of lo!i2 GeV/c 2 ; 
under the assumption that the 'e error' results from uncalculated higher-order pQCD 
contributions the effects represented by these variations are already 'counted', and we 
have not added them in quadrature with the other theoretical uncertainties. Their 
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inclusion does not change the central rrib(Mz) value and increases the total theoretical 
uncertainty to t?;fg GeV/c 2 . 



Our result is in agreement with that from ||13|| . The latter measurement has a 
significantly smaller theoretical uncertainty. For the Durham algorithm alone we would 
obtain an uncertainty of similar size, but our study of six different jet algorithms has 
revealed additional systematic effects which warrant further investigation. 
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Figure 1: The R\/ R% ds ratios measured by SLD for each of the six jet finding al- 
gorithms (horizontal bands) compared with the predicted dependence on the running 
6-quark mass, m b (M z )- 
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